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ABSTRACT 

We report a dynamical measurement of the mass of the brown dwarf GJ 802B using aperture- 
masking interferometry and astrometry. In addition, we report the discovery that GJ 802A is itself 
a close spectroscopic non-eclipsing binary with a 19 hour period. We find the mass of GJ 802 B to 
be 0.063 ± 0.005M Q . GJ 802 has kinematics inconsistent with a young star and more consistent with 
the thick disk population, implying a system age of ^10 GYr. However, model evolutionary tracks for 
GJ 802B predict system ages of ~2 GYr, suggesting that brown dwarf evolutionary models may be 
underestimating luminosity for old brown dwarfs. 
Subject headings: stars: low-mass, brown dwarfs 
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1. INTRODUCTION 

The boundary between stars and brown dwarfs is de- 
fined as the mass at which the luminosity of old objects 
is just dominated by hydrogen burning. This boundary 
is predicted to occur at a mass between 0. 07-0.072 M P) at 
[Fe/H1=0 and ~0.092Mq at [Fe/H]=-3 (jChabrier et all 
I2M01: l&ilrrows et al.ll2T)0lh . 

However, this boundary, and the theoretical relation- 
ships that predict effective temperature and luminosity 
as a function of mass and age, are largely untested by 
observations. The observations that are required to test 
these models are dynamical mass measurements of bi- 
nary and multiple-star systems, combined with accurate 
photometry and distance determinations, preferably at 
known age. Many stars betwee n 0.1 and 0.2 Mp, now 
have accurate dynamical masses (jSegransan et al.l l2000). 
but objects with masses between the hydrogen burning 
limit and 0.1 M© so far do not have accurate (< 10%) 
mass and luminosity determinations. 

Dependence on theoretical models has led to contro- 
versy about the mass of brown dwarfs, particularly when 
objects are plausibly near th e canonical planetary -mass 
boundary of ~13Mj up (e.g. iLuhman et al.l [2007) . In 
the few cases where the dynamical mass of a brown 
dwarf has been me asured (|Zanatero Osorio et al.l 120041 ; 
IStassun et all l2006f ), precision is either inadequate to 
truly constrain models, or the model fits have to 
make unusual assum ptions such as non-coeval systems 
(IStassun et al.ll2007h . 

GJ 802 is a M5 field dwarf system at ~16pc that was 
discovere d to have a brown dw arf component through as- 
trometry (|Pravdo et all 2005) . The subsequent detection 



of the compa nion, GJ 802B, w ith aperture-masking in- 
terferometry (jLlovd et al.ll2006T ) made this an ideal target 
for dynamical mass determination, due to the system's 
~3 year period. The high contrast (~100:1) and small 
separation (< 100 mas) of GJ 802AB make this system 
too difficult for direct imaging observations, but is in an 
ideal range for interferometric detection. 

This paper reports aperture-masking interferometry 
detections of GJ 802B at six epochs in three colors, 
enough to make good measurements of the orbit, and 
the discovery that GJ 802A is itself a close spectroscopic 
M dwarf binary. Section [2] describes the astrometric, 
interferometric, spectroscopic and photometric observa- 
tions that go into the orbit determinations. Section [3] 
describes the orbit of the close binary GJ 802Aab and 
Section U describes the constraint on the dynamical mass 
of GJ 802B from the orbit of GJ 802 AB. In Section [5] we 
compare our results to theoretical predictions and Sec- 
tion [5] has a summary and discussion. 

2. OBSERVATIONS 

The observations on which this paper is based consist 
of seeing-limited astrometry, aperture-masking interfer- 
ometry, infrared and visible spectroscopy and photome- 
try. The astrometry observations, in addition to parallax 
and proper motion, provide the orbit of the GJ 802Aab 
pair about the GJ 802AB center of mass. The aperture- 
masking interferometry constrains the wide binary orbit 
by resolving GJ 802B. The spectroscopy, in the infrared 
and visible, provides the spectroscopic GJ 802 Aab orbit, 
and the photometry is used to search for eclipses of this 
inner pair. Each observation set will be discussed in turn. 
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TABLE 1 

ASTROMETRY OBSERVATION SUMMARY FOR GJ 802 



Date ( U 1 ) 


JU-2450000 


AHA (mas) 


ADec (mas) 


1998-07-06 


1000.75 


0.00 


0.00 


1999-07-20 


1379.75 


901.50 


1806.99 


1999-09-23 


1444.75 


1002.89 


2092.99 


2000-07-13 


1738.75 


1780.05 


3510.93 


2001-07-02 


2092.75 


2631.23 


5150.19 


2002-07-07 


2462.75 


3513.06 


6908.15 


2003-09-09 


2891.75 


4496.65 


8932.57 


2004-07-12 


3198.75 


5276.92 


10367.24 


2004-09-15 


3263.75 


5371.75 


10654.40 



2.1. Astrometry Observations 

The astro metry observa t ions i n this paper come di- 
rectly from iPravdo et alj (|2005l ). These observations 
were part of the Stellar Planet Survey (STEPS) program, 
with observations made using a custom visible camera 
mounted at the unfolded Cassegrain focus of the Palo- 
mar 200 inch telescope. As the raw astrometry was not 
published in that paper and have not been made avail- 
able to us, we extracted the astrometric data from the 
figures in the paper, and then added the motion due to 
the proper motion and parallax given in the paper's ta- 
ble. These extracted values are presented in Table [TJ 
and include all ep ochs where there are at least 2 mea- 
surements given in IPravdo et all (|200l . We give all of 
these epochs equal weight. The error estimated from the 
scatter within individual epochs is 1.6 mas, however we 
assume errors of 1.7 mas, to give a final fitted x 2 of 1.0 
(see Section 0}. 

2.2. Aperture-masking inter 'ferometry 

The technique of non-redundant aperture masking has 
been well-established as a means of a chieving the full 
diffraction limit of a single telescope (e.g . lMichelsonlll920t 
iBaldwin et al]|1986l : [Tuthill et al.ll2000D . It involves plac- 
ing a mask with an array of holes of non-redundant spac- 
ing in the pupil-plane of a telescope, and analyzing the 
recorded images as interference fringes on a number of 
discrete baselines. Recently, our group has used this 
technique behind adaptive optics at the Palomar 200 inch 
and Keck telescopes to increase the magnitude limit of 
the technique beyond that achievable in a seeing-limitcd 
speckle mode. The reason for the technique's success 
over direct adaptive-optics imaging is that the calibra- 
tion is independent of structure of the wavefront over 
scales larger than a single sub-aperture, but it still pre- 
serves the angular resolution of the full aperture. 

Aperture-masking observations of GJ 802 were made 
using bot h the PHARO camera of the Palomar 200 inch 
telescope (|Havward et al.ll2001[ ) and the NIRC2 camera 
of the Keck II telescope. The Palomar apert ure- masking, 
mode of the PHARO camera is described in iLlovd et alJ 
(2006). The NIRC2 aperture- masking mode is similar, 
although the masks are placed in the filter wheels of the 
NIRC2 camera rather than a designated pupil wheel as 
is the case with PHARO. The masks create fringe pat- 
terns in the image-plane of each detector. Our primary 
observable, closure-phase, is extracted from these fringe 
patterns and models are fit to the closure phases. A sum- 
mary of all aperture-masking observations, along with 
the detected binary properties, are given in Table [2] 



For the NIRC2 and PHARO aperture-masking exper- 
iments, both 9 and 18 hole masks are available. The 9 
hole mask results in 36 simultaneous baselines and 28 in- 
dependent closure-phases. The 18 hole mask results in 
153 simultaneous baselines and 136 independent closure- 
phases, but has half the throughput of the 9 hole mask 
and spreads the light over 4 times as many pixels as the 
9 hole mask. 

All observations were done with a 9 hole mask in H 
and K bands, except the June 2007 J CO nt observation for 
which a 18 hole mask was used. In the absence of sky 
rotation, the field-of-view of the aperture-masking ex- 
periment is given by X/2Bs, where B$ is the shortest 
baseline in the mask. Outside this field-of-view, the po- 
sition angle and separation become ambiguous without 
multiple exposures at different sky rotations. Although 
the 9 hole mask was well suited to H and K observations 
at both Keck and Palomar, the ^73 mas separation of GJ 
802B would have been outside the nominal 9 hole mask 
field of view in J band. 

In a similar manner, the aperture masking inner work- 
ing angle is given by A/2B^, where is the longest 
baseline (near to the full aperture size). At this separa- 
tion, the maximum closure phase signal is equal to the 
contrast ratio of the binary (i.e. 0.01 radians or 0.6 de- 
grees for a 100:1 binary). In the case of the PHARO 
observations of June 2004 and October 2006, these nom- 
inal inner detection limits are 43 and 57 mas respectively. 
Near this inner edge, the contrast ratio and separation 
for model-fitting are degenerate. Despite their large er- 
rors on separation, these detections are therefore reliable. 
We explicitly list the correlations between separation and 
contrast in Table 2. 

We have made several improvements to the analysis 
pipeline and observi ng procedur e since the original de- 
tection published in ILlovd et al.l (|2006h . The most im- 
portant improvement has been in the observations them- 
selves: we have generally been much more careful in us- 
ing calibrators of similar brightness, colors and position 
in the sky. This procedure was not used carefully for the 
discovery epoch of September 2004, and we find that the 
astrometry depends on which of several unsuitable cal- 
ibrators are used in the analysis. Therefore, we choose 
not to use this epoch in our analysis here, and it is not 
presented in Tabled 

We have windowed the data prior to Fourier- 
transforming with a tighter function than previous pub- 
lications to minimize residual chip-based effects (e.g. of 
bad pixels). The window size is of the form exp(— r 4 ) 
with a full-width half-maximum of 1.6X/Dh, where Dh 
is our aperture-mask hole diameter projected onto the 
primary mirror. This tighter window also has the effect 
of spatially-filtering the interferometric data. The raw 
frames from the PHARO camera are taken in a mode 
where all (non-destructive) reads are saved. This allows 
us to split the data in post-processing into a variety of 
integration times, where smaller integration times have 
less atmospheric and AO system noise but more read- 
out noise. We now carefully choose an optimal number 
of sub-reads to analyze data of a given brightness: for 
GJ 802 we split the data into sub-frames spaced by 2 
reads, which have integration times of 862 ms each. 
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TABLE 2 

Aperture-masking observation summary for GJ 802 



Date (UT) 


JD-2450000 


Instrument 


Filter 


Separation 


Position Angle 


Contrast a 


Sep. /Contrast 








(mas) 


(degrees) 


(B/A) 


correlation 


2004-06-06 


53163.0 


PHARO 


H 


55. 1± 21.4 


16.4± 6.8 


0.013±0.006 


0.87 


2006-06-23 


53910.0 


NIRC2 


Ks 


85.8± 3.5 


200. 9± 2.2 


0.015±0.003 


0.09 


2006-10-10 


54018.7 


PHARO 


H 


94.6± 9.7 


207.6± 7.4 


0.009±0.003 


-0.13 


2006-10-10 


54018.7 


PHARO 


Ks 


81. 1± 24.5 


196. 7± 4.9 


0.011±0.003 


0.75 


2007-05-31 


54252.0 


PHARO 


CH4S 


77.4± 10.3 


24.1± 5.2 


0.009±0.002 


-0.26 


2007-06-05 


54257.0 


NIRC2 


Jcont 


73.3± 2.9 


21. 5± 2.1 


0.008±0.002 


0.04 


2007-06-06 


54258.0 


NIRC2 


Kp 


72.6± 1.7 


19.2± 1.2 


O.OlOiO.001 


0.01 


2007-07-31 


54312.9 


PHARO 


H 


86.2± 5.8 


24.0± 3.5 


0.010±0.002 


-0.40 


2007-08-29 


54341.8 


PHARO 


H 


85.4± 6.0 


20.2± 3.2 


0.011±0.002 


-0.44 
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2.3. Infrared Spectroscopy 

The NIRC2 camera at the Keck II telescope was used 
on UT 2006-08-10 in a spectroscopic mode for the pur- 
pose of attempting to detect the spectrum of GJ 802B 
directly. The observations were made in the H-band with 
a bandpass from 1.52 to 1.62 microns. We used a grism 
with a resolving power of 17580 per pixel and a 4 pixel 
wide slit, giving a spectral resolution of ~4000. These 
observations did not succeed in their primary purpose, 
mainly because poor seeing forced us to use H instead of 
J-band, where the spectral differences between a mid-L 
and mid-M dwarf were not large enough to give a clearly 
detectable GJ 802 B signal in the data. 

It was noticed, however, that lines in the spectra of 
GJ 802 were double. To calibrate the spectra, we first 
used Xe and Ar lamps to calibrate the wavelength scale, 
which was fitted with a second-order polynomial. We cal- 
culated the model atmospheric and grism transmission 
by taking spectra of the F8V star HD 136118 and divid- 
ing b y a template spect rum of HR 4375 (a GOV star) 
from iMever et all (j 1998) . In reducing these spectra to 
their two components, GJ 777B was used as the spectral- 
type standard and the radial velocity standard, as its 
wide companion GJ 77 7A has a precise radial velocity of 
-45.350±0.004kms- 1 (|Naef et al.1 12001 . and should be 
within lkms -1 from GJ 777B's radial velocity due to 
potential orbital motion. As a check of our calibration 
procedure, we used Xe and Ar lamp calibration of the 
NIRC2 grism wavelength scale to measure the radial ve- 
locity of GJ 777B based on two Al I and one Ca I lines 
to be —41 ± 5kms _1 , consistent with GJ 777A's radial 
velocity. 

Least square fits were made to the continuum sub- 
tracted spectra of GJ 802, based on a model made up of 
the sum of two shifted GJ 777B spectra. We chose this 
technique rather than a cross-correlation because the two 
spectra were not separated well enough to give clearly 
separate peaks in the cross-correlation. The errors in the 
velocity difference come straight from the least-squares 
fitting process, where the uncertainty in the GJ 802 spec- 
trum was set to 1.3% of the mean flux so that the model 
fit had a reduced \ 2 of 1. An uncertainty on the veloc- 
ity difference based on deviations from a linear fit to the 
measured velocities would be 2.5 km s" 1 . The results of 
this fitting process arc given in Tabled The magnitude 
difference in _ff-band derived from the fit is zero within 
errors, which suggests that the components have nearly 
equal masses. 

It is difficult to estimate the uncertainties on the abso- 
lute velocity calibration (not the primary purpose of our 
observations), so we assign an error of 5kms _1 based 
on the absolute calibration of GJ 777B's radial veloc- 
ity against the Xe and Ar lamps. There is good reason 
to expect that the absolute calibration of these observa- 
tions is significantly worse than the relative calibration 
between the Aa and Ab spectra, because the point spread 
function delivered by the AO system was not guaranteed 
to be well-centered on the slit. However, as both stars 
had the same point spread function, with their maximum 
possible spatial separation of 0.7 mas corresponding to 
~lkms _1 , the relative calibration between Aa and Ab 
spectra should be good compared with the 2kms _1 sta- 
tistical error in the velocity difference. 



TABLE 3 

Infrared Radial Velocities for GJ 802Aab 



JD- 2450000 


RV A6 


RV A a 


ARV 




3956.9219 


-6.9±5.0 


-93.5±5.0 


-86.6±2.0 


0.15±0.06 


3956.9336 


-13.2±5.0 


-86.1±5.0 


-72.9±2.0 


0.19±0.06 


3956.9375 


-12.9±5.0 


-83.7±5.0 


-70.9±2.0 


0.24±0.06 


3956.9453 


-16.3±5.0 


-83.6±5.0 


-67.3±2.0 


-0.03±0.07 



a Difference in apparent H-band magnitudes as derived from the fit of 
two shifted GJ 777B spectra. See text. 



2.4. High- Resolution Optical Spectroscopy 

We obtained multi-epoch spectra of GJ 802 using the 
East-Arm Echelle (EAE) on the Palomar Hale 5m tele- 
scope. The EAE is a high-resolution spectrograph, capa- 
ble of achieving spectral resolution of R ~ 30000. Most of 
our observations span a wavelength range of 4000-10000 
A, but the observations taken in December 2006 (before 
official commissioning) span a significantly bluer range 
(3500-8000 A). Most of our spectra were binned in the 
spectral direction in order to reduce read noise, so the 
effective resolution for our observations was R ^20000. 
Since GJ 802 is very red, none of our spectra have any 
useable signal short-ward of ^5500 A. 

We observed GJ 802 at four epochs on 2006 December 
15 and six more epochs spanning 2007 May 06-08. We 
also observed the late- type stars GJ 51 (in December) 
and GJ 581, GJ 686, and GJ 699 (in May) as spectral 
type standards. All observations of GJ 802 used integra- 
tion times of 300-900s, while our standard star observa- 
tions used a wider range of integration times since some 
were significantly brighter. 

We bias-subtracted, flat-fielded, and extracted our 
spectra using standard IRAF tasks. Wavelength cali- 
bration was achieved with respect to a Thorium-Argon 
lamp that was observed at the beginning of each night; 
preliminary tests suggest that telescope flexure results in 
wavelength calibration variations of no more than ^0.015 
A (~0.5 km s _1 ) in the vicinity of the Ha emission line at 
6563 A. In Figure [TJ we show a representative segments 
of several spectra for GJ 802 around the Ha wavelength 
range. These spectra demonstrate the double-lined na- 
ture of GJ 802; this plot also shows that the Ha emission 
line strengths are not constant, but vary on a time-scale 
of days. 

All of our science spectra exhibit Ha emission, so we 
have directly determined the component radial veloci- 
ties from the line centroids. For epochs where the lines 
were not clearly resolved, we fit the spectra with a pair 
of Gaussian emission lines with the same FWHM as the 
resolved measurements. We list all of these radial veloc- 
ities, including heliocentric corrections, in Table [5J we 
also list our observed radial velocity standards in Ta- 
ble \5\ The mean radial velocity of the standards was 
3.0 km s" 1 higher than the radial velocities derived from 
the Thorium- Argon lamp calibration, with a dispersion 
of 1.7kms -1 (standard error on the mean 0.5kms _1 ). 
We have subtracted this offset from all radial velocities 
reported in Table E] and Table |U 

We did not use any cross-correlation techniques to de- 
termine overall RV fits because we are still attempting 
to characterize the wavelength and flux calibrations of 
the instrument, including limits on potential variability; 
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TABLE 4 

High-Resolution Spectroscopy of GJ 802. 



JD - 2450000 


RV A6 


RVa 


a 


A RV 


4084.57847 


-107.7± 


5.0 


10.2± 


5.0 


-117.9± 


2.8 


4084.59097 


-111.8± 


5.0 


15. 7± 


5.0 


-127.5± 


2.8 


4084.60347 


-118.2± 


5.0 


18. 0± 


5.0 


-136. 2± 


2.8 


4084.70625 


-116.9+ 


6.0 


19. 4± 


6.0 


-136. 3± 


4.2 


4084.71389 


-118.7± 


6.0 


18. 5± 


6.0 


-137.2± 


4.2 


4084.72014 


-113.2+ 


7.0 


16. 6± 


7.0 


-129. 8+ 


7.1 


4226.94971 


-102.3+ 


2.0 


23. 2± 


2.0 


-125. 5± 


2.8 


4227.90039 


-99.3+ 


2.0 


17.9± 


2.0 


-117.2± 


2.8 


4228.84473 


-21.4± 


2.0 


-59. 8± 


2.0 


38.4+ 


2.8 


4228.90332 


9.7+ 


2.0 


-86. 7± 


2.0 


96.4+ 


2.8 



TABLE 5 

High-Resolution Spectroscopy of Standards. 



Star 



JD - 2450000 



RV 



Literature RV 



GJ 


51 


4084.66134 


-10.9±3 


-7.7 a 


G.I 


393 


4226.70127 


+4.6±3 


+8.3 b 


G.I 


476 


4226.75556 


+29.6±3 


+31. 6 a 


G.I 


686 


4226.83443 


-11.4±3 


-9.5 b 


G.I 


476 


4227.73785 


+30.6±3 


+31. 6 a 


G.I 


514 


4227.74175 


+11.2±3 


14.6 b 


G.I 


526 


4227.74530 


+13.0±3 


16.0 a 


G.I 


581 


4227.74907 


-13.8±3 


-9.4 b 


G.I 


393 


4228.63097 


+4.9±3 


+8.3 b 


G.I 


387.2 


4228.64279 


-16.8±3 


-16.5 a 


G.I 


699 


4228.82439 


-116.9±3 


-110.5 b 



1 IGizis et all d2002fl b INidever et~aT1 J2002D 





655B B5B0 6562 6561 

Wavelength (ang) 

Fig. 1. — High-resolution spectra for GJ 802 centered on the Ha 
emission line. Both components show resolved Ha emission, but 
the Ha emission line strengths appear to be variable on a time-scale 
of days. 



given the long time baseline of our dataset, any improve- 
ment in the precision would yield only very minor im- 
provements in the orbital ephemerides. 



2.5. Photometry 

Given that the GJ 802Aab binary had such a short pe- 
riod, and preliminary analysis of the orbit suggested that 
inclination was possibly high, we used the Palomar 60- 
inch robotic telescope to search for eclipses. The Palomar 
60-inch telescope has a single instrument, a CCD imager 
with an 11x11 arc- minute field-of view and a selection 
of broad and narrow-band filters. We took series of im- 
ages in g, i and H-a filters around the predicted times 
of eclipses based on the radial-velocity orbit (when the 
velocity difference was zero). The H-a filter was used to 
search for an eclipse of the chromospheric emission. 

The observations in g and i filters spanned 30 min- 
utes on either side of a predicted eclipse at Julian day 
2454311.889, and the observations in H-a spanned 30 
minutes on either side of a predicted eclipse at Julian 
day 2454315.870. These eclipse times were based on a 
preliminary orbital fit (see Section [3]). As the radial ve- 
locity amplitudes of Aa and Ab were roughly equal (see 
the velocities in Table[4]), we can assume that the masses 
are roughly equal and that the primary and secondary 
eclipses would have roughly equal depths. 

We computed aperture photometry for GJ 802 on the 
calibrated images, using the median photometric varia- 
tion of 10 nearby field stars to compute the sky trans- 
mission in each frame. No eclipses were found in any 
filters. For the g and i filters, the RMS photometry scat- 
ter was 0.01 and 0.003 mag respectively with sampling at 
1.5 minute intervals. For the H-a filter, the scatter was 
0.01 mag, and the sampling at 1 minute intervals. The 
excess in our filter due to the H-a emission was 22.7% of 
the continuum, calculated by comparison with photom- 
etry from an off-line narrow-band filter and consistent 
with our spectroscopic observations. From these data we 
can place 2a upper limits for a continuum eclipse depth 
of 0.005 mag and for a H-a eclipse depth of 0.01 mag. 

3. ORBIT OF THE CLOSE PAIR GJ 802A 

We chose to fit the radial velocities of GJ 802Aab with 
a circular orbit only, because the 19 hour period of the 
system is 10 times shorter than the canonical cutoff of 7 
days f or tidal circularization in low-mass main-sequence 
stars (Za hn fc Bouchetl[l989l ). Figure [2] shows the best 
fit equal-mass circular orbit with the measured radial 
velocities, and Figure [3] shows the best fit circular orbit 
with the measured radial velocity differences. We made 
this fit by first examining by eye the 2007 May radial ve- 
locities, concluding that the orbital semi amplitude was 
between 135 and 155 kms" 1 , with a period of approxi- 
mately 0.8 days. We then fit to the absolute value of the 
velocity difference by using a grid search. We searched 
the range of 0.7 to 0.9 days with 10 -5 days spacing for pe- 
riod, the full range of epoch (modulo half a period) with 
10~ 3 days spacing, and total velocity amplitudes of the 
two components Kl + K2, between 135 and 155 kms -1 
with 4kms~ 1 spacing. The reduced \ 2 f° r the final fit 
is 0.70 with 11 degrees of freedom. No other (aliased) 
minima have reduced x 2 values less than 6. We assigned 
signs to the radial velocity differences only after finding 
this fit, as GJ 802 Aab is so close to equal brightness that 
it was difficult to tell which spectrum was GJ 802Aa and 
which was GJ 802Ab. 

Best fit orbital parameters are given in Table [5] It is 
not clear from these observations which component is the 
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TABLE 6 
Radial Velocity Solution for 
GJ 802Aab 



Parameter 


Value 


Epoch (HJD) 


2454140.530±0.001 


Period (d) 


0.795340±0.000003 


Kl + K2 


149.1±1.5 kms" 1 


Kl 


73.8±1.4 kms" 1 


K2 


75.4±1.4 kms" 1 


(M A a+M A b)sin 3 i 


0.273±0.008 M Q 


a sin(i) 


2.343±0.024 R Q 


Inclination 


77<i<83 



50 




- 1 fin I 

0.0 0.2 0.4 0.6 0.8 1.0 

Phase 

Fig. 2. — Best fit orbit for GJ 802Aab, showing the Aa velocities 
(solid line, diamonds) and the Ab velocities (dashed line, triangles). 



200 




-200 1 -\ 

CO 0.2 0.4 0.6 0.8 1.0 

Phase 

Fig. 3. — Best fit orbit for GJ 802Aab, showing the difference in 
velocities between the two components. 

more massive, as the mass ratio Mj^ > /M j ^ a (= K1/K2) is 
0.98±0.03. This ambiguity is seen in Figure^ The mass 
ratio and center of mass radial velocity are constrained 
best by the May 2007 observations, which are best cali- 
brated and where the two spectra swap their positions, 
and the absolute calibration of the radial velocities is 
most certain. These data have significantly smaller er- 
rors in Figure [2] and are at phase ~0.03, 0.1, 0.65 and 
0.85. Including a ~3 kms -1 uncertainty for the orbital 
motion of GJ 802Aab with respect to the GJ 802AB cen- 
ter of mass, we obtain a radial velocity for the system of 
-42±4kms _1 . 

The inclination limit in Table [5] comes from the fol- 
lowing analysis of both the lack of eclipses and a model- 



dependent mass estimate for the system. Based on a 
stellar radius of_0.16i?©, predicted by the models of 
iBaraffe et al.l ()1998f ) for components of mass ~0.14M© 
(applicable to an edge-on system with equal-mass com- 
ponents and (M J 4 a +M J 4b)sin 3 i from Table [6]), the eclipse 
would have lasted a maximum of 50 minutes. For a graz- 
ing eclipse with 5.8% of one star eclipsed, the eclipse 
would have lasted 25 minutes: well within our 60 minute 
observing window (Section l2.5p . A grazing eclipse occurs 
at a projected separation of twice the stellar radius, or 
~0.32 Rq. As no eclipse occulting even 1% of the surface 
of one of the stars was found in any filter (Section 12. 5| , 
the inclination of GJ 802Aab is limited to < 83 degrees. 
Given that (M A a+MAb)sin 3 i = 0.273±0.008, this in turn 
limits the total mass of GJ 802A to > 0.279 ± 0.008. 

The models of Baraffc et al.l ()1998l ) have been verified 
to predict the K-band mass-luminosity relationship for 
field dwarfs correctly to withi n ~5% in mass for d warfs 
of mass greater than 0.1 M© (|Derfosse et al"1l2004[ ). For 
stars of mass ~0.11-0.18M©, this relationship is inde- 
pendent of age at the 1% level for ages between 0.5 and 
10 Gyr. This age range is applicable to the likely GJ 802 
age of >6 Gyr, based on its activity (jPravdo et al.ll2005f ) 
and kinematics (Section[5]). This relationship is also rela- 
tively steep: at a mass of 0.13 MQ, a 6% error in K-band 
flux from the distance of iPravdo et alJ (|2005l ) translates 
to only a 2.5% error in mass. Therefore, the uncer- 
tainty in the mass is likely dominated by the K-band 
mass-luminosity relationship (at least until the models 
are tested at higher precision) , and we assign a standard 
deviation of 5% to the model-predicted mass. These 
models predict masses for GJ 802Aa a nd GJ 802Ab 
of 0T 34Mq, using a parallax of 63 mas IPravdo et al.l 
(2005) and 2MASS photometry. In turn, this places a 
maximum-mass constraint on GJ 802Aab at ~ 2<r of 
0.295 M©, limiting the inclination of GJ 802Aab to be- 
tween 77 and 83 degrees. Note that the if -band lumi- 
nosity - mass relationship is little affected by metallic- 
ity, changing by < 1% if we were to use the models of 
iMontalban et all (l2000h at [Fe/H]=-1. 

4. DYNAMICAL MASS OF GJ 802B 

The full solution for the GJ 802AB orbit requires 16 
parameters: 5 for parallax and proper motion, 7 for the 
orbital solution of B with respect to Aab, 1 for the photo- 
center semi-major axis, and 3 for the contrast of B with 
respect to Aab in J, H and K bands. The reason that 
contrast had to be added in to the orbital solution was 
that many of the parameters in Table [5] had strong de- 
generacies. To fully explore this large parameters space 
around the best-fit solution, we used a Ma rkov Chain 
Monte Carlo technique (e.g. iBremaudl 1 1 999h . a method 
that has often been used in a stronomy for cosm ologi- 
cal parameter estimation (e.g. iKnox et al.ll2001h . Key 
advantages of this technique are the ability to easily in- 
clude the covariance matrix of the data (in our case the 
aperture-masking fits), and the ability to easily calculate 
the posterior probability function of derived parameters, 
such as the mass of GJ 802B. 

The probability that a particular set of parameters 
is contained in the final chain is proportional to the 
likelihood function, which is proportional to exp(x 2 /2). 
Due to the correlations between derived parameters from 
aperture-masking, x 2 is n °t just the sum of normalized 
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deviates, but makes use of the covariance matrix of the 
data: 

X 2 = (m(0)-d)*C- 1 (m(0)-d). (1) 

Here d is the vector of data values, m is the model for 
these data based on parameters 9, C is the data covari- 
ance matrix and * represents a transpose. Covariances 
between different epochs were assumed to be zero: only 
the derived separation, position angle and contrast for 
a single epoch of aperture-masking data had non-zero 
covariances. 

We used a Markov Chain of length 1.4 x 10 5 , with a 
1.4 x 10 4 burn- in time. The best orbital solution had a re- 
duced x 2 of 1.0 because we chose the error in the STEPS 
data to be 1.7mas in each axis. As in iPravdo et al.l 
( 2005) , we used a 2 ± 1 mas conversion from relative to 
absolute parallax. 

We first conducted an exploratory unconstrained fit 
to the data, but found that the errors on parallax and 
orbital semi-major axis were too large to give a useful 
total mass. This fit with single parameter ler confidence 
limits is given in the first column of Table [7J and has a 
mass for GJ 802Aab that is too high by ~l.lcr from the 
mass derived from the absolute magnitude of GJ 802Aab. 
Therefore, we fixed fixed the masses of GJ 802Aa and 
GJ 802 Ab to be 0.134 Mq, ba sed on their absolut e K- 
magnitudes and the models of iBaraffe et al.l (|1998[ ) (see 
Section[3]) . Note that the same value of parallax of 63 mas 
was used in both Section [3] and here in calculating the K- 
band absolute magnitude. This fit is given in the central 
column of Table [0 

Finally, we note that we can not really "fix" the mass 
of GJ 802A as it has an uncertainty. We therefore cal- 
culated an apriori constraint on the mass of GJ 802Aab 
using the following: the mass-K magnitude relationship 
with its assumed 5% RMS error (see Section[3]); the value 
of Msin 3 (i) from Table El and an orbital orientation as- 
sumed random with i < 83 degrees. Together, these 
constraints give a mass of 0.2775 ± 0.0082 for the mass 
of GJ 802Aab. The assumed likelihood curve resulting 
from these two constraints is given in Figured) it can be 
seen that a Gaussian likelihood curve is a good approxi- 
mation for our prior knowledge of Ma ■ The knowledge of 
this mass was added to the Markov-Chain Monte Carlo 
as an apriori constraint to give our preferred values in 
the final column of Table ® The final GJ 802B mass of 
0.063±0.005M Q is the most accurate mid-late L-dwarf 
mass reported in the literature so far. 

Although the Markov-Chain Monte Carlo technique is 
an excellent technique for accurately calculating poste- 
rior probability distributions of derived parameters, it 
can be difficult to intuitively understand the magnitude 
of derived uncertainties from the chain output. There- 
fore, we will examine the origin of the 8% uncertainty on 
the G J 802B mass individually. The mass of the B com- 
ponent directly relates to the mass of the A components 
by the ratio of astrometric and orb ital semi-major axes 
a and a (e.g lMcCarthv et al.lll988h : 

M B = -^—M A . (2) 

1 — a/a 

Therefore, a 4% uncertainty on the mass Ma, a 6% 
uncertainty on a and a 2% uncertainty on a combine 




0.22 0.24 0.26 0.28 0.30 0.32 
M (Solar Masses) 

Fig. 4. — The apriori constraint on the mass of GJ 802A, from the 
K-band luminosity (dotted line), the radial velocity orbit (dashed 
line), and the joint likelihood (solid line). 

to give a 8% uncertainty on Mb- Although in the con- 
strained fit, the value and uncertainty for a is influenced 
by the Ma constraint, this argument for the value and 
uncertainty of Mg changes little if the values for a and 
a from the unconstrained fit are used. 

The apparent photometry corresponding to this pre- 
ferred fit is given in Tabled] The apparent photometry 
for the s ystem's combined light comes from Simbad and 
2MASS (|Cutri et al.l l2003f) . The derived absolute pho- 
tometry corresponding to individual components of the 
system (apportioning equal flux to the two components 
Aa and Ab) is also given in this table. This photometry 
corresponds to a sp ectral-type of ^L5-L7 for GJ 802B 
(jKnapp et al.ll200l . 

The orbital fit for GJ 802AB is shown in Figure 02 
with the photo-center astrometry measurements scaled 
by the ratio of total mass to GJ 802B mass, and with 
parallax and proper motion removed so that all points 
can be plotted at the same scale. This plot also aids in 
developing an intuitive feel for our quoted mass error. 
The error in the ratio between the mass of GJ 802B and 
the total mass is the ratio of the photo-center and orbital 
semi-major axes (similar to Equation [5]): which we claim 
has a 7% error. 

5. COMPARISON WITH MODELS 

As discussed in (Pra vdo et al.ll2005T ). the activity of 
GJ 802 (e.g. Figure[T]), imply that the system is >6Gyr 
old. The kinematics of GJ 802 also imply that the system 
is old. The total proper motion of nearly 2 arcseconds 
per year is 30 times the parallax, meaning that the tan- 
gential velocity of GJ 802 is ~30AU per year, which is 
~140kms~ 1 . Based on the astrometric fit reported in 
Table [7J and the mean radial velocity from Tables [3] and 
IH we have calculated the space motion of GJ 802 using 
the gal_uvw routine from the IDL Astronomy User's Li- 
br ary 1 . This routine is in t urn ba sed on the mathematics 
in Uohnson fc Soderbloml (|1987t) . The (U, V, W) space 
motion calculated to be (134, — 51, 21)kms~ 1 , which is 

1 http://idlastro.gsfc.nasa.gov/ 
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TABLE 7 

ASTROMETRIC SOLUTION FOR GJ 802AB 



Parameter 



Unconstrained Fit 



Fixed 



Constrained M A 



Proper Motion (RA) 
Proper Motion (Dec) 
Parallax 

Orbital semi-major axis (a) 

Photocenter semi-major axis (a) 

Epoch (JD) 

Period 

Eccentricity 

Argument of periapse 

Longitude of ascending node 

Inclination 

Total Mass 

Mass of GJ 802Aab 

Mass of GJ 802B 



877.7±1.0 mas yr" 1 

1722.11±0.34 mas yr" 1 

62.2±1.8 mas 

97.5±5.0 mas 

17.4±1.2 mas 

2453040±14 

1105± 9 d 

0.40±0.08 

90.4± 4.3 

21. 8± 1.4 

83.7± 3.0 

0.426±0.078M o 

0.351±0.066M o 

O.O76±O.O13M 



877.7±1.0 mas yr" 1 
1722.2±1.0 mas yr' 1 
63.9±1.3 mas 
92.4±2.0 mas 
17.2±0.9 mas 
2453028±15 
1105± 9 d 
0.35±0.05 
89. 0± 5.5 
22. 7± 1.4 
81. 1± 3.1 
O.329±O.OO4M 
0.268 a 

O.O61±O.OO4M 



877.7±1.0 mas yr -1 

1722.1±1.0 mas yr" 1 

63.5±1.3 mas 

92.9±2.1 mas 

17.2±1.0 mas 

2453031±13 

1104± 9 d 

0.35±0.05 

89. 6± 4.8 

22. 0± 1.3 

82. 7± 3.0 

O.343±O.O12M 

O.28O b ±O.OlOM 

O.O63±O.OO5M 



a Fixed from a mass - K magnitude relationship (see text). Mass with error-bars included as apriori information 
in the Monte-Carlo fits. 



TABLE 8 
Absolute Photometry for GJ 802 



Band 


m A 


m B 


M Aa (=M Ab ) 


M s 


ya 


14.67 




14.47 




J 


9.57±0.02 


14.75±0.27 


9.34±0.05 


13.74±0.28 


H 


9.07±0.02 


14.13±0.09 


8.83±0.05 


13.14±0.10 


K 


8.76±0.01 


13.61±0.08 


8.53±0.05 


12.62±0.08 



a GJ 802B is assumed to have negligible contribution to V-band. 

(134, —39, 28) km s _1 with respect to the local standard 
of rest. The lag in V and very high U space velocity 
is quite inconsistent with the thin-disk population. For 
example, the HIPPARCOS sample of nearby early M 
stars has a U dispersion of 32km/s, and the sample ear- 
lier than spectral- type F5, representative of a ~1 GYr 
old population , has a U dispersion of only 22kms _1 
(|Mignardll2000fl . The 39 km s" 1 V lag of GJ 802 is unusu- 
ally small for a Popu lation II halo star, but is consistent 
with the thick disk (jCasertano et al.lll990l ). GJ 802 is 
very unlikely to be a runaway star that is confused with 
the thick disk population, because any dynamical inter- 
action capable of giving it a ~ 100km s^ 1 peculiar veloc- 
ity would also break the wide binary GJ 802 AB apart. 
Therefore, these kinematic properties place GJ 802 at an 
age of ~10 GYr, an d almost certainly older than 3 GYr 
(jBensbv et al.ll2003l) . 

We can compare these ages to the ages derived from 
the modelling of GJ 802B. Figure [5] shows the absolute 
magn itudes of GJ 802B co mpared to the DUSTY mod- 
els of Baraffc et al.l (|2002t ). The colors of these models 
are clearly too red, a well known property of these mod- 
els for objects that are of mid-L or later spectral type 
(jBaraffe et al.ll2003[) . The predicted age is between 1 
and 5 GYr in age in all bands. We chose not to plot 
the results from the COND models from the same group 
(|Baraffe et alj r2003) because these models are much too 
blue and really only applicable to T dwarfs. 

Although these models clearly do not represent 
GJ 802B, mostly due to difficulties in opacity treatment 
and dust in brown dwarf atmospheres, it is expected 
that luminosity and effective temperature are a much 
more certain prediction of brown dwarf evolutionary 
models than fluxes in specific bands (|Burrows fc Liebertl 
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Fig. 5. — Best fit orbit for GJ 802 AB (green), with aperture 
masking points in blue and scaled seeing-limited astromctry in red 
(see text). 



11993 iChabrier et alJl2000h . Abetter comparison for the 
GJ 80 2B photometry is th erefore possible with the mod- 
els of lBurrows et al.l (|2006j ) that correctly model the basic 
colors of L dwarfs and the L to T transition. 

At log(<?)=5.2 (appropriate for a ~0.06 brown 
dwarf), the best fitting model from lBurrows et all (|2006D 
has absolute magnitudes Mj = 14.05, Mjj = 13.10 and 
M K = 12.78. This model has T cS = 1400 K with a to- 
tal luminosity of 2.9 x 10~ 5 L Q . This model Luminosity 
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Fig. 6. — Absolute magnitudes of GJ 802B compared to the DUSTY models at 0.5, 1 and 5 GYr. The comparison is for J-band (left 
panel), H-band (center panel) and K-band (right panel). Model ages are 0.5 GYr (dotted line), 1 GYr (solid line) and 5 GYr (dashed line). 



and temperature can now be compared with evolutionary 
models, pl acing a 0.063 Mq GJ 802B at 1.9 GYr accord- 
ing to the iBurrows et~aT1 (|1997l ) track. GJ 802B would 
only be consistent with an old (i.e. >5GYr) object if 
its mass were 0.07 Mq, which is inconsistent with our 
dynamical mass at 1.5cr. 

6. CONCLUSIONS AND DISCUSSION 

GJ 802 is a triple system with component masses of 
~0.14, 0.14 and 0.063 Mq. The inner pair of equal-mass 
stars, GJ 802Aab, has an orbital period of 0.795 days 
and an inclination between 74 and 83 degrees. The outer 
component, GJ 802B, has a ~3 year orbit and a mass of 
0.063 ± 0.005 M Q . The inclination of the GJ 802B orbit, 
^83 degrees, is consistent with co-planarity. The most 
promising way to decrease the mass uncertainty would 
be to obtain an accurate radial velocity curve of the 3 
year orbit, with an expected amplitude of ^3km/s. 

The brown dwarf "desert" is generally defined as a 
lack of low-mass companions around solar-type stars, 
but a similar lack of very unequal mass-ratio companions 
has b een found around very low-mass stars (jClose et al.1 
2003). However, as a triple system, GJ 802 no longer 
fits in this category. In fact, there are a large number 
of brown dwarfs known in triple or higher order multiple 
systems. This entirel y consistent with dynamical star- 
formation simulations (|Delgado-Donate et alj|2004j ). but 
may also relate to the details of the fragmentation pro- 
cess as a cloud with relatively high angular momentum 
collapses. Measuring the mutual inclination (i.e. degree 
of co-planarity) of this and other low-mass triples with 
long-baseline interferometry could help to determine the 
mechanism responsible for producing triples like GJ 802. 

The kinematics of GJ 802B place it in the thick disk, 
at an age of ^10 GYr, while the model-derived age for 
GJ 802B is ^2 GYr. Although this discrepancy is only 
significant at the 93% confidence level so far (1.5c), we 
will list several possibilities for reconciling the discrep- 
ancy: 1) GJ 802 Aab could actually have a total mass 
of >0.30Mq, consistent with the astrometry, but be 
underluminous by ~ 30% in K-band compared to field 
stars. GJ 802B would then be placed right on the sub- 



stellar boundary at ~0. 07 Mq, a naively unlikely posi- 
tion given the very small sample of astro metric STEPS 
binar ies from which GJ 802 was taken (|Pravdo et al.l 
2005). 2) GJ 802Aab could have an apparent total mass 
of >0.30Mq because of an additional low-mass compo- 
nent in a ~ 30 day, just stable orbit. This could also 
increase our dynamical mass of GJ 802B to ~0. 07 Mq. 
3) GJ 802 could actually be <2 GYr old but have experi- 
enced a unique dynamical past that gave it a high space 
velocity without tearing the wide binary apart. Or, 4) 
models for old brown dwarfs are systematically under- 
predicting luminosities. These model errors could relate 
to, e.g. the effects of magnetic fields that ca n hinder heat 
flow in brown dwarfs (jChabrier et al.ll2007t ). Given that 
GJ 802B is currently the only ~Gyr or older brown dwarf 
with an accurate dynamical mass, we find this fourth 
possibility most likely: a hypothesis that will be testable 
within the next few years as more field brown dwarfs 
have accurate mass determinations. 
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